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Abstract: In this paper we propose a straightforward method for understanding 
the thermodynamics of black holes in de Sitter space, one that will allow us to study 
these black holes in a way that is analogous to the anti-de Sitter case. As per usual, we 
formulate separate thermodynamic hrst laws for each horizon present in the spacetime, 
and study their thermodynamics as if they were independent systems characterized by 
their own temperature. That these systems are not entirely independent and various 
thermodynamic quantities in them are in fact ‘correlated’ is reflected by the fact that 
their thermodynamics can be captured by a single Gibbs free energy-like thermody¬ 
namic potential. This quantity contains information about possible phase transitions 
in the system and allows us to uncover a rich phase structure for de Sitter black holes. 
In particular, we discover reentrant phase transitions for Kerr-dS black holes in six 
dimensions, a phenomenon recently observed for their six dimensional AdS cousins. 
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1 Introduction 

Since the early 1970’s it has been known that the characteristic properties of black 
holes such as the surface gravity, horizon area, mass, etc. are related by a law anal¬ 
ogous to the first law of thermodynamics [1, 2]. Ever since then the subject of black 
hole thermodynamics has been an active area of research that has been both exciting 
and insightful. In particular, significant attention has been devoted to the study of 
asymptotically anti-de Sitter (AdS) black holes. For these black holes, thermodynamic 
equilibrium is straightforward to define and a wealth of interesting phenomena emerge. 
For example, a thermal radiation/large black hole first order phase transition was found 
by Hawking and Page for the Schwarzschild-AdS black hole [3] and the gravitational 
analogue of the liquid/gas phase transition was observed in the phase diagram of the 
charged AdS black hole [4-7]. More exotic AdS black hole spacetimes may even demon¬ 
strate reentrant phase transitions [8] or admit a tricritical point in their phase diagram 
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[9] (see [10-12] for a review of these developments). The study of phase transitions 
in AdS black hole spacetimes is especially interesting due to their dual conformal held 
theory (CFT) interpretation through the AdS/CFT correspondence [13]. 

More recently, a holographic duality between gravity in de Sitter space and con¬ 
formal held theory has been proposed [14], spurring interest in the thermodynamics of 
black holes in asymptotically de Sitter (dS) spacetimes. Asymptotically dS black holes 
are also of direct interest in cosmology. Thermodynamics in dS space is, however, bur¬ 
dened with issues not present in the AdS case. The existence of a cosmological horizon 
in addition to the black hole horizon means that the system associated with an observer 
living between the horizons is in a non-equilibrium state—such an observer would hnd 
herself in a thermodynamic system characterized by two temperatures. Additionally, 
the lack of a global timelike Killing vector outside the black hole prevents one from 
dehning a good notion of the black hole mass (see [15, 16] for recent developments in 
this direction). 

One way to deal with the thermodynamics of asymptotically dS black hole space- 
times is to formulate several separate thermodynamic hrst laws, one for each “physical” 
horizon present in the spacetime [17-20] (see also [21]). Specihcally, let us consider a 
general rotating dS black hole with several U{1) charges in d dimensions. Typically, 
such a black hole admits three horizons, located at real positive radii r which are de¬ 
termined from the horizon condition, say /(r) = 0. The cosmological horizon (we use 
subscript c) is located at the largest positive root Tc, the black hole horizon (denoted by 
subscript h) corresponds to the second largest positive root r^, and the inner horizon 
(with subscript i) corresponds to the third largest positive r* (if it exists). The hrst 
laws for the three horizons then read [20] d 

m = mst + - tLW + Vt5P, ( 1 . 1 ) 

k j 

m = -TJS, + - aDsr + - <s>i) 6 <^ + v,sp , ( 1 . 2 ) 

k j 

SM = -T.SS. + J])!!? - - t’JiQ’ + V^SP . (1.3) 

k j 

Here, M represents a quantity that would be the ADM mass in asymptotically AdS and 
hat cases and for this reason we simply refer to it as “mass”.^ The horizon temperatures 

^ The thermodynamics of the inner black hole horizon might seem rather formal. However, only 
after it is included in the thermodynamic description, thermodynamic quantities demonstrate a certain 
type of ‘universality’, which is speculated to provide some information about the microscopies of black 
holes, see e.g. [22-24]. 

^In the dS case such a quantity is “conserved in space” (rather than in time) due to the spacelike 
character of the Killing field dt in the region near infinity [25]. 
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Tb, Tc and Ti are all defined to be positive. Sb, Sc and Si denote the horizon entropies, 
ri’s and J’s denote the angular velocities and momenta respectively, <h’s and Q’s stand 
for the electric potentials and charges, and the pressure P is related to the positive 
cosmological constant A according to 


^ _ (d- l)(d-2) 

Stt IQtiP 


(1.4) 


as it would be for a perfect fluid stress-energy tensor. Treating the cosmological con¬ 
stant as a thermodynamic pressure in this way is an idea that has recently received 
much attention [12]. I 4 , I 4 and Vi are the thermodynamic volumes—the quantities 
thermodynamically conjugate to P —given by 


' \dP) 




Vb={ 




V 






(1.5) 


Starting from (1.1)-(1.3), it is possible to formulate the so called ‘subtracted’ hrst 
laws. For example, for an observer in between the cosmological and the black hole 
horizon, the difference of (1.1) and (1.2) gives 


0 = Tb5Sb + Tc5Sc + - V6P , (1.6) 

* j 


where V stands for the subtracted volume of the ‘observable universe’, V = Vc — Vb > 0. 
Although such subtracted laws provide useful information, in this paper we directly 
work with relations (1.1)-(1.3). 

The above three laws are accompanied by the corresponding Smarr-Gibbs-Duhem 
formulae 

= ns, + ^ - ^V,P, (1.7) 

j k 

= -TA + ^ -aDn- , ( 1 , 8 ) 

j k 

= -T.S. + ^ + ^(nf - Sli) J" - , (1.9) 

j k 


which can be derived from the corresponding hrst laws (1.1)-(1.3) via the dimensional 
scaling argument, see e.g. [26]. 

Note that the three hrst laws (1.1)-(1.3) are not truly independent and the various 
thermodynamic quantities in them are in fact ‘correlated’ and determined from the 
same number of independent parameters as in the AdS case. For example, consider the 
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Schwarzschild-dS black hole spacetime in four dimensions, with the line element given 
by 


= —fdt^ + 


T 


+ r^dQ 


2 ! 


^ ^ 2M ^ SvrPr^ 


( 1 . 10 ) 


and dfl^ denoting the element on the unit 2-sphere. By hxing the cosmological constant 
P and the mass M, the cosmological and black hole horizon radii {k = b, c) are given 
by the largest (uc) and second largest (r;,) positive roots of /(n^) = 0. (Note that 
there is no inner horizon in this case.) Having hxed these radii, both temperatures 
are then uniquely determined by = |/'(rfc)|/(47r), as are the entropies Sk = vrr^, 
and thermodynamic volumes 14 = In other words, the complete thermodynamic 

description of all horizons is hxed by choosing M and P, exactly as in the Schwarzschild- 
AdS case. This immediately generalizes to arbitrary dS black holes. As we shall, see 
this will allow us to construct a Gibbs free energy-like quantity that captures behavior 
of all three dS horizons and is in some sense a continuation of the AdS Gibbs to negative 
pressures and temperatures. 

The aim of this paper is to study the thermodynamic behavior of asymptotically dS 
black holes and, in particular, to investigate if such black holes admit phase transitions 
similar to their AdS cousins [3-12]. 

One approach to such a task is that of the effective temperature [27-32].^ In this 
approach one concentrates on an observer who is located in an ‘observable part of the 
universe’, in between the black hole horizon and the cosmological horizon. For such an 
observer a new effective hrst law is imposed such that the system is assigned a ‘total 
entropy ’ that equals the sum of the black hole horizon and the cosmological horizon 
entropies, see e.g. [35, 36], 

S = St + S,, ( 1 . 11 ) 


while the parameter M is interpreted as internal energy, and the parameter V = 14 — 14 
as the observable volume. This then dehnes the effective temperature and pressure— 
quantities thermodynamically conjugate to S and V. Interestingly, this prescription 
leads to a Schwarzschild-dS black hole analogue [37] of the famous Hawking-Page 
transition [3]. However, in general the physical meaning or the positive-dehniteness 
of the effective temperature and pressure are not obvious. Moreover, applying the 
same prescription to more complicated dS black hole spacetimes results in complicated 
phase diagrams whose interpretation is not clear [37]. 

In this paper we propose an alternative path towards the thermodynamics of dS 
black holes. Namely, we do not attempt to assign the system a single effective temper¬ 
ature. Instead, we follow the following steps. 

•^See also [34] and [33] for recent alternative approaches. 
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1. We study the thermodynamics of all three physical horizons—treating them as 
‘independent thermodynamic systems ’ in equilibrium, characterized by their own 
temperature, thermodynamic behavior, and hrst law (l.l)-(1.3)d 

2 . Since the thermodynamic quantities on these horizons are (at least for a given 
classical solution) correlated, one can assign to the whole system a ‘single’ Gibbs 
free energy-like quantity G. This quantity, which is in some sense a continuation 
of the Gibbs free energy of AdS black holes to negative pressures and negative 
temperatures, captures the behavior of all three horizons.® 

3. We assume that we have a mechanism (external force) at our disposal that allows 
us to tune thermodynamic parameters of the horizons.® 

4. Possible instabilities/phase transitions at each horizon, encoded in quantity G, 
are taken as a sign of instability/phase transition of the whole dS system.^ 

Effectively, our procedure boils down to Ending the Gibbs free energy G. By study¬ 
ing its properties one can understand the thermodynamic behavior of each horizon as if 
it were an independent thermodynamic system. Putting all the information together, 

“^Let us stress that we do not claim that the entire dS black hole system is in thermodynamic 
equilibrium. That this is not the case is, for example, reflected by the fact that the corresponding 
Euclidean instanton is not regular unless the ratio of horizon temperatures is a rational number. Since 
each horizon is in general at a different temperature, there will be a heat flow in the system that might 
eventually, if the system were not disturbed, equilibrate the horizon temperatures. However, in what 
follows we assume that such a process would happen at a time scale which is much longer (due to the 
smallness of h) than a typical timescale of thermodynamic phase transitions induced by an ‘external 
force’ considered below. 

^The definition of G is not an additional assumption but rather a mathematical convenience that 
directly follows from 1) and the correlation of thermodynamic quantities at different horizons. In our 
study we discard the possibility of thermal phase transitions governed by Euclidean instantons that 
would ‘change the character’ of the classical gravitational solution under consideration, e.g. [38]. 

®This is a standard assumption made in black hole thermodynamics in order to describe possible 
phase transitions. For example, in the famous Hawking-Page transition [3] it is assumed that an 
infinite reservoir can ‘pump’ energy into the system, changing radiation to a large black hole phase 
as the temperature of the system increases. A novel feature that seems to emerge in dS black hole 
spacetimes is that a change in thermodynamic parameters of one horizon necessary leads to a change 
of thermodynamic parameters of the other two horizons. Although strange at a first sight, exactly 
the same correlation silently takes place in asymptotically flat or AdS black hole spacetimes where 
the inner horizon is assumed to adjust accordingly to the change of the black hole horizon. For a 
discussion of this phenomenon for the charged AdS black holes see App. A. 

^That phase transitions can happen for systems with temperature gradients and far from thermal 
equilibrium is well documented by our everyday experience. Consider, for example, a typical ice fishing 
story: while we boil water to make tea and get warmer, our ice fishing hole freezes over. 
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one can then infer possible black hole phase transitions. Treated in this way we discover 
that dS black holes admit phase transitions and complicated phase diagrams, similar 
to the AdS case. For example, as we shall see, one hnds a reentrant phase transition 
analogous to [8]. 

The rest of the paper is organized as follows. In the next section we pedagogically 
illustrate the above procedure using the simple example of four-dimensional charged 
dS black holes, studying the thermodynamics of each horizon separately and showing 
how this can be described by a single thermodynamic function G. Sec. 3 is devoted 
to rotating dS black holes in four and higher dimensions; in particular the reentrant 
phase transitions are discovered for the 6-dimensional case within a certain range of 
parameters. Sec. 4 is devoted to conclusions. App. A describes a similar construction 
applicable to AdS black holes. App. B contains technical results regarding the rotating 
dS black holes in all dimensions. 


2 Thermodynamics of three horizons: charged dS case 

In this section we show how the thermodynamics of the three de Sitter horizons can be 
captured by a single Gibbs free energy, and how such a quantity encodes information 
about possible phase transitions of the de Sitter black hole system. For concreteness 
and technical simplicity this is illustrated on a ‘test spacetime’ of the charged dS black 
hole in four dimensions. 


2.1 Charged dS black holes 


The charged dS black hole line element takes the following form: 


ds'^ = —fdt^ H —— 


r'^dQl , 


2M g^p^2 

J f 9 o 


( 2 , 1 ) 


where M stands for the mass and Q for the electric charge. The Schwarzschild-dS 
metric is recovered by setting Q = 0. The cosmological r^, black hole r^, and inner p 
(for g 7^ 0) horizons are located at Vc > ri, > Vi > 0, given by f{rk) = 0. The cases 
when two or more horizons coincide are special and are partly discussed below. The 
thermodynamic quantities associated with the three horizons are given by, see e.g. [20], 


n = 

1 

dTrrft 

-F SvrPrft - 

Q\ 

rir 

Sb = Tirl , 


g 

n ’ 

Vb 

4 3 

= -Tirf , 

3 ^ ’ 


Tc = 

1 

f 1 -I- STiPrj — 


Sc = vrr^, 

<hc = 

g 

Vc 

4 3 

= -vrr^, 



AtITc 


^ ’ 



rc ’ 

3 ’ 


T- = 
1 

( 

dvruj 

^1-1- SuPrl — 


Si = Txr] , 


Q 

Ti 

Vi-- 

4 3 

= -Tirf , 

3 * ’ 

(2.2) 
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and obey the three hrst laws (1.1)-(1.3), as well as the corresponding Smarr relations. 

Several special cases of charged dS spacetimes are of importance. The extremal 
black hole limit corresponds to the case where the inner horizon coincides with the black 
hole horizon, rt = r^. It can be easily shown that in this case we have Tb = Ti = 0. 
The second important case is the Nariai limit in which = r^, and consequently, 
Tb = Tc = 0. Finally, in the limit r;, —)■ 0 we recover the pure dS case, while in the limit 
Tc —)■ oo the asymptotically flat black hole solution is obtained. We refer to [38-42] for 
more details on these special cases. 

2.2 Black hole horizon 

Let us start our discussion by studying the behavior of the black hole horizon. The 
corresponding thermodynamic system is characterized by the black hole horizon tem¬ 
perature Tb, while the hrst law (1.1) now reads 

5M = Tb6Sb + <I>b5Q + Vb6P. (2.3) 

The structure of this formula suggests that the mass M should be treated as the 
gravitational enthalpy, M = Hb, similar to the AdS case [26]. In fact, the above hrst 
law is exactly identical to the hrst law of AdS black holes, see e.g. [10], with the 
exception that P = — A/Stt is now negative. For this reason one might expect that 
the thermodynamic behavior of the dS black hole horizon can be obtained by simply 
‘continuing’ the AdS thermodynamics to negative pressures. The quantity analogous 
to the Gibbs free energy G = G{Tb, P, Q) is given by the Legendre transformation of 
the enthalpy M, 

G = Hb-TbSb = M-TbSb. (2.4) 

Note that Tb is manifestly non-negative and vanishes in both the Nariai limit and the 
extremal black hole limits. 

Having dehned the Gibbs free energy, we can display its behavior as a function 
of the temperature Tb and reveal the associated thermodynamic behavior. For our 
test spacetime this is displayed by thick black curves in Fig. 1. In the charged case we 
observe two branches of black holes (meeting at a cusp), compared to one branch in the 
Q = 0 limit. This behavior of G can be understood as the continuation of the known 
behavior of the Gibbs free energy in the AdS case, displayed in the hgure by dashed 
red curves, and, apart from the ‘Nariai cut-off’, it is in many respects similar to the 
behavior of the Gibbs free energy of asymptotically flat spacetimes, c.f. dashed black 
curves. Whereas the charged AdS case (see App. A for more details) admits a critical 
point characterized by a critical pressure Pc = l/QGvrQ^, and below such pressure one 
observes a small black hole/large black hole phase transition, in the case of charged dS 
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Figure 1. Gibbs free energy: black hole horizon. The Schwarzschild (left) and the 
charged with Q = \ (right) cases are displayed for various pressures. The thick black curve 
corresponds to the dS case with P = —0.0045, the dashed black curve to the asymptotically 
flat case with P = 0, and the red dashed lines to the AdS case with P = 0.0015 (below 
Pc) and P = 0.0040 (above Pc), respectively. Clearly the dS case can be understood as a 
continuation of the AdS case to the region of negative pressures. Since the critical point (in 
the charged case) occurs for positive pressures, there is no criticality in the dS case. Note 
however that the upper branch in between the Nariai limit and the cusp is thermodynamically 
unstable as it has negative specific heat and higher Gibbs free energy than the lower branch 
that terminates at the extremal limit. 


black holes no such behavior exists. There is no swallow tail and the behavior of G 
is relatively simple: two branches (the upper one being thermodynamically unstable) 
meet at a cusp. 

To probe the local thermodynamic stability, one can study the specific heat. For 
the black hole horizon this is naturally defined as 


a 


dH, 


m 


P,Q 


dM 


m 


P,Q ' 


(2.5) 


As expected, for the Schwarzschild-dS case Cb is negative and vanishes in the Nariai 
and dS limits, see left Fig. 2. In the charged dS case, the black hole branch with lower 
Gibbs (in between the extremal black hole and the cusp) has positive Cb and is locally 
thermodynamically stable, at the cusp Cb suffers from an infinite jump, and is negative 
for the upper branch of black holes (in between the cusp and the Nariai limit). In the 
Nariai and extremal limits Cb vanishes, see left Fig. 3. 













Figure 2. Specific heats: Schwarzschild-dS case. We display the two specific heats: 
Ch (left) and Cc (right) for P = —0.009 as a function of x = rhjrc. Cb is non-positive and 
vanishes in both the Nariai (x = 1) and dS (x = 0) limits. Cc is non-negative and vanishes 
in the Nariai limit. 





Figure 3. Specific heats: charged dS case. We display the three specific heats: Cb (left), 
Cc (middle) and Ci (right) for Q = 1 and P = —0.009. Cb suffers from an infinite jump at x 
corresponding to the cusp (C) of the Gibbs free energy and vanishes in both the Nariai (N) 
and extremal (E) limits. Cc vanishes at the Nariai limit and otherwise is manifestly positive, 
Ci vanishes in the extremal limit and otherwise is negative. 

2.3 Cosmological and inner horizons 

Let us next investigate the thermodynamic behavior of the cosmological horizon, char¬ 
acterized by the temperature Tc- Staring at Eq. (1.2), we realize that the TcdSc term 
has an opposite sign to what is customary in order to treat M as the enthalpy. To fix 
this, we set He = —M < 0, and write (note the flipped sign of the work terms) 

6Hc = Tc6Sc - ^eSQ - Vc6P . ( 2 . 6 ) 
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Figure 4. Gibbs free energy: cosmological horizon. The Schwarzschild-dS (left) and 
the charged dS with Q = 1 (right) cases are displayed for P = —0.0045. In both instances Gc 
monotonically decreases and does not indicate the presence of any thermodynamically unsta¬ 
ble regions or phase transitions. However, for the charged case (right) a thermodynamically 
unstable region, indicated by the dashed black curve, exists for Tc G (0, Tq). This region 
corresponds to the upper thermodynamically unstable branch of G for the black hole horizon 
displayed in Fig. 1; Tq is the temperature of the cosmological horizon that corresponds to the 
cusp of the black hole horizon Gibbs free energy. 


In other words, we identify the following Gibbs free energy: 

G, = G,{T,, P, Q) = -M- T,Sc . (2.7) 

For our test spacetimes, this is displayed in Fig. 4. On its own, the behavior of Gc 
seems ‘boring’ and does not indicate the presence of any thermodynamically unstable 
regions or phase transitions. However, this is misleading as, similar to the AdS case 
discussed in App. A, one has to ‘coordinate’ the behavior of the cosmological horizon 
with the behavior of the black hole horizon, i.e., to consider G and Gc simultaneously. 
Consequently, for the charged dS black hole a part of the Gc graph becomes thermody¬ 
namically unstable (displayed by dashed black curve), as it corresponds to the upper 
(thermodynamically unstable) branch of G, c.f. Figs. 1 and 4. 

Finally, for the charged dS case the thermodynamic behavior of the inner horizon 
should also be taken into account. This is formally similar to that of the cosmological 
horizon, with replaced by r*, resulting in the definition Gi = Gi(Ti, P,Q) = —M — 
TiSi, and is displayed in Fig. 5. 

One can also study the specihc heats of the cosmological and inner horizons. These 
are now defined as 


dHc 


dM 


^ m 


dM 

dTc 

P,Q 

dTc 

P,Q ’ 

* m 

P,Q 

dTi 
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Figure 5. Gibbs free energy: inner horizon. The Gi — Ti diagram is displayed for the 
charged dS black hole with Q = 1 and P = —0.0045. 

and displayed for our test spacetimes in Figs. 2 and 3. We find that Cc is manifestly 
positive (apart from the Nariai limit where it vanishes) for both charged and uncharged 
dS black holes. G* is only dehned for the charged dS case; it is negative (apart from 
the extremal limit where it vanishes) indicating the local thermodynamic instability of 
the inner horizon. 

2.4 Unified description of all horizons 

To treat all the horizons simultaneously, we note that 


Gcivc, P, Q) = -G{rb -)■ Tc, P, Q ), 

Tc{rc, P, Q) = -Thin Tc, P, Q), (2.9) 


and similarly for the inner horizon. This means that, similar to the AdS case discussed 
in App. A, the information about the thermodynamics of all horizons is in fact encoded 
in G, dehned by Eq. (2.4), provided we extend its validity to “a// admissible radiP ri, 
and correspondingly to negative temperatures T^. Such negative temperatures then 
correspond to the positive temperatures of the cosmological and inner horizons and in 
these regions the thermodynamic equilibrium occurs for the global maximum (rather 
than minimum) of G. 

Specihcally, let us consider the charged dS black hole with Q 7 ^ 0. The thermody¬ 
namics of all three horizons is encoded in G = G(T, P, Q), given by 


G 


^ 2 3 ^ 3g2 

- ttGP H- 

4 3 4 r 


T 


--h 2rP — 

47rr 


47rr^ 


( 2 . 10 ) 
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Figure 6. Gibbs free energy: all horizons. The Gibbs free energy G is displayed 
for all possible horizons: inner (blue), cosmological (red), and black hole (black) for the 
Schwarzschild-dS (left) and the charged dS with Q = 1 (right) cases for P = —0.0045; dashed 
curves correspond to regions that are thermodynamically unstable. Note that the positive 
x-axis corresponds to T = Tf, and the negative x-axis to T = —Tc or T = —Tj temperatures. 
For a given state of the black hole this allows one to read temperatures of all horizons. 
For example, the cusp occurs for the black hole horizon temperature ~ 0.016; 

Tj = Tig ~ 0.083 and Tc = Teg ~ 0.013 are the corresponding temperatures of the inner and 
cosmological horizons. 

Let us fix P and Q and denote by the corresponding radius of the extremal black 
hole and by the radius of the Nariai black hole. We further denote by the minimal 
radius of the inner horizon, which occurs for the Nariai limit, and by tm the maximal 
radius of the cosmological horizon, which occurs for the extremal black hole. We then 
parametrically plot G = G(T) for the range r G (rm,rM), being careful to interpret 
the thermodynamic quantities for different ranges of r. For r G {rm,rE), G = —Gi and 
T = —Ti represent inner horizon quantities and the global maximum of G corresponds 
to a thermodynamically preferred state. For r E {rE,rN), G = Gb and T = Tb represent 
black hole horizon quantities while G is minimized by the thermodynamically preferred 
state. For r G {rE^rM), G = —Gc and T = —Tc represent cosmological horizon 
quantities and the preferred state again corresponds to the maximum of G. This is 
illustrated in right Fig. 6. 

The case with Q = 0 is similar, with the exception that there is no longer an inner 
horizon. Consequently, the range for r is r G (0,rM) and we have r G (0,r7v) for the 
black hole horizon quantities and r G for the cosmological horizon quantities, 

as displayed on the left of Fig. 6. 

Having dehned the Gibbs free energy G that describes the thermodynamic behavior 
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of all involved horizons, we can identify all possible thermodynamically unstable regions 
and phase transitions present in the system. For example, consider the charged dS black 
hole with Q = 1 and P = —0.009 whose thermodynamics is displayed in Fig. 7. On the 
left side we depict the behavior of the three horizon temperatures, displayed as function 
of a; = Tb/rc- Note that only values of x between extremal (E) and Nariai (N) bounds 
are admissible, x G {xm = Oii fhe right side the Gibbs free energy is displayed. 

Let us hrst concentrate on the black hole horizon. Choosing some admissible black hole 
horizon radius = Tb^ (corresponding to x = Xq) in the lower black hole branch, we 
can ask what would happen, for example, with the dS system provided we increased 
the black hole size Vb a little bit. From the graphs one can immediately read all the 
temperatures as well as how they are going to change as we increase r^, following the 
small black arrows. Note also that as long as changes in are small (until we meet the 
cusp) we move along a lower black hole branch that has lower Gibbs free energy than the 
upper branch and hence is thermodynamically preferred. However, crossing the cusp 
would move the system in the upper black hole branch which is thermodynamically 
unstable. Consequently, we observe unstable regions (dashed) for not only the black 
hole horizon but also the cosmological (red) and inner (blue) horizon branches. This is 
how the information about the thermodynamics of each horizon interplays and builds 
towards the understanding of the thermodynamic behavior of the whole system. 

3 Rotating dS black holes 

In the previous subsection we have studied the thermodynamics of charged dS and 
Schwarzschild-dS black holes. We observed that both examples are relatively ‘boring’. 
Whereas the Gibbs free energy for the Schwarzschild-dS system is completely smooth, 
the Gibbs free energy of charged dS black hole at least admits a cusp and conse¬ 
quently thermodynamically unstable regions, however no interesting phase transitions 
are present in either case. Interesting phase transitions, indicated for example by the 
presence of a swallow tail, might exist for more complicated dS black hole spacetimes, 
possibly in higher dimensions. In this section we uncover such behavior by considering 
rotating dS black holes in four and higher dimensions. 

The general rotating dS metrics in all dimensions and the corresponding thermody¬ 
namic charges are gathered in App. B. Similar to the charged case, we have up to three 
horizons and the corresponding Nariai and extremal black hole limits. To investigate 
the thermodynamic behavior we study the Gibbs free energy G = G{T, P, P) in the 
canonical (hxed P) ensemble. The corresponding expression is slightly different for 
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Figure 7. Thermodynamics of charged dS black holes. Left: The horizon temperatures 
are displayed as functions of x = r^frc- We observe 3 important values of x: the extremal 
limit (E), the Nariai limit (N), and a cusp (C) for which the black hole horizon temperature 
reaches its maximum, xq is chosen as an illustration point; it determines Ti,q,Tcq and 
for the r.h.s. graph. The arrows illustrate what happens with horizon temperatures as the 
size of the black hole rb increases from rf,p. Namely, the black hole horizon and inner horizon 
temperatures locally increase whereas the cosmological horizon temperature decreases. Right: 
The corresponding Gibbs free energy, displaying the thermodynamic behavior of all three 
horizons. In both graphs P = —0.009 and Q = 1. 


odd and even dimensions, see App. B. For example, in even dimensions we have 


G = 


^d-2 

Stt 


N 

n 

k=l 


+ a| rl - 


^ 1 
1 

1=1 


+ Oi 


1 + r^jl 

2r 


2 -, 


T 


r(l — r^/P) 

2^r 


N 

E 

2=1 


+ a] 


l + ryP 

dvrr 


(3.1) 


where N = 2N + 2, and a* are understood as functions of the angular momenta, 
Qi = ai{Jj) through formulae (B.4) and (B.2). 


3.1 Behavior in four and five dimensions 

It turns out that (similar to the AdS case) the behavior of the Gibbs free energy of 
4d and 5d rotating dS black holes is qualitatively similar to the behavior of G of the 
four-dimensional charged-dS black hole. Namely, we observe no interesting features 
apart from a cusp of divergent specihc heat. This is displayed in Fig. 8. In the five¬ 
dimensional case we have checked that the same qualitative behavior is obtained for 
various ratios q = J 1 /J 2 of the two angular momenta. As we shall see in the next 
subsection, a new feature appears in six dimensions. 
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Figure 8. Thermodynamics of rotating dS black holes in fonr and five dimensions. 

Left figure displays the G — T diagram for the four-dimensional Kerr-dS black hole for P = 
—0.0025 and J = 1. Right figure displays the G — T diagram for a doubly-spinning Kerr-dS 
black hole in five dimensions for P = —0.05 and q = J 1 /J 2 = 0.5. In both instances we 
observe qualitatively the same behavior as in the 4-dimensional charged-dS case. 

3.2 Six dimensions: reentrant phase transitions 

In the six-dimensional case the behavior of G becomes more interesting. Namely, in 
Fig. 9 we display the Gibbs free energy for the doubly spinning Kerr-dS black hole in 
d = 6 dimensions for fixed angular momentum ratio q = J 1 /J 2 = 0.112.® The com¬ 
plicated behavior of G (in the black hole horizon region) indicates the presence of a 
reentrant phase transition, exactly analogous to the one observed for AdS black holes 
in [8]. Specifically, as the black hole horizon temperature Tf, monotonically increases, 
we observe a reentrance of the small black hole phase as follows. Below = Ti a, 
(lower left) branch of small black holes is preferred, at Ti a first order phase transition 
occurs and the (lower middle) large black hole phase becomes temporarily dominant. 
This continues until at Ti, = Tq the large black hole branch terminates. Increasing 
the temperature even further results in a zeroth order phase transition as the systems 
‘jumps back’ to the (lower right) branch of small black holes. One can easily check that 
in this case the branches with lower Gibbs (for a given temperature) are locally ther¬ 
modynamically stable as they have positive specific heat. In other words, we observe a 
small/large/small black hole reentrant phase transition, similar to the AdS case.® 

®As with the AdS case, the particular behavior of G depends crucially on the choice of ratio Ji/J 2 . 
For example, for a given q = 0.112, we can find the reentrant phase transition in the range of pressures 
P G (0,Pc) where Pc ~ —0.0022. 

®Note, however, that in the AdS case the reentrant phase transition occurs for singly spinning 
Kerr-AdS black holes in six dimensions, whereas the doubly spinning ones may feature a tricritical 
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Figure 9. Reentrant phase transition. The black hole horizon Gibbs free energy is 
displayed for the doubly spinning Kerr-dS black hole in d = 6 dimensions for hxed angu¬ 
lar momentum ratio q = J 1 /J 2 = 0.112 and P = —0.001. We observe a typical behavior 
characteristic of the reentrant phase transition, previously observed for the AdS black holes 
[9]. Namely, as the temperature increases the system undergoes small/large/small black hole 
reenrant phase transition: as the temperature increases, the first order phase transition from 
small to large black hole at T = Ti is followed by a zeroth order phase transition at T = Tq 
as the system ‘jumps back’ to the small black hole branch. 

4 Conclusions 

The subject of black hole thermodynamics continues to be one of great importance in 
gravitational physics. Over the past hfty years many of the studies in this held have 
concentrated on the thermodynamics of asymptotically hat and asymptotically anti 
de Sitter black holes, the latter emerging in light of the AdS/CFT correspondence. 
However, the thermodynamics of asymptotically de Sitter black holes possesses new 
difficulties and remains relatively unexplored. One of the primary reasons is that de 
Sitter black hole spacetimes are essentially non-equilibrium multi-temperature systems; 
each horizon is characterized by its own temperature and a hrst law-like relation. Al¬ 
though such laws are not independent and various quantities in them are correlated, 
it is by no means obvious, despite several attempts in the literature, that a reasonable 
equilibrium thermodynamic description based on the ehective temperature approach 
can be consistently formulated. This however does not mean that one has to give up 

point and small/intermediate/large black hole phase transition [9]. 
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on the thermodynamics of these systems. 

In this paper we proposed an alternative strategy for stndying the thermodynamics 
of de Sitter black holes. Namely, we formulated several thermodynamic hrst laws (one 
for each horizon, as usual) and studied their thermodynamics as if they were indepen¬ 
dent thermodynamic systems characterized by their own temperature. The internal 
correlation of these systems then enabled us to dehne a single Gibbs free energy-like 
quantity G that captures their individual thermodynamic behavior, as well as contain¬ 
ing information about possible phase transitions of the whole de Sitter system once 
properly combined. In this way we were able to uncover the rich phase structure of de 
Sitter black holes. Namely, although we found that the thermodynamic behavior of the 
four-dimensional Schwarzschild-dS and charged dS black holes is relatively uninterest¬ 
ing, we discovered the presence of reentrant phase transitions for Kerr-dS black holes 
in six dimensions. 

As a simple consistency check we also dehned and investigated the horizon spe- 
cihc heats. In particular, we found that for the studied cases whenever a possibility of 
multiple black hole branches emerges (for a given temperature) the branch with pre¬ 
ferred Gibbs has also positive specihc heat and is locally thermodynamically stable. Is 
our alternative proposal for studying the thermodynamics of dS black holes internally 
self-consistent? And if so, how many of the phase transitions and critical phenomena 
observed for the AdS black holes hnd their analogue in the dS case? These remain 
interesting open questions left for future investigations. 
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A AdS analogue of liquid/gas phase transition: inner horizon 

In this appendix we re-investigate the thermodynamic analogue of a liquid/gas phase 
transition present in the charged AdS black hole spacetime [4-7]. The purpose is 
twofold, i) To show that a “multi-horizon description” similar to the one studied in 
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Figure 10. AdS analogue of liquid/gas phase transition. Left: The Gibbs free en¬ 
ergy G is displayed for the charged AdS black hole for Q = 1 and P = 0.003; we observe 
the characteristic swallow tail behavior: the thermodynamically unstable configuration of the 
outer black hole horizon is displayed by the dashed red curve, temperature Th = Tq corre¬ 
sponds to the first-order phase transition. Right: The corresponding P — Tf, phase diagram 
is reminiscent of the liquid/gas phase transition. 

the main text is already (silently) present in the AdS case, ii) To address the following 
question overlooked in the existing literature: What happens to the inner horizon when 
the outer horizon experiences a Van der Waals-like phase transition? 

To this end, we consider a charged asymptotically AdS black hole in four dimen¬ 
sions, given by the metric 

ds'^ =-fdt"^+ ^ Pr'^dnl, / = !- —+ ^ + ^, (A.l) 

and identify the negative cosmological constant with positive pressure, according to 
P = [d — l){d — 2)/ (Ibvr/^). The black hole horizon is located at the largest > 0 for 
which f{rf,) = 0. The thermodynamic quantities then read 

Tb = (l + SvrPrfc - , St = Tirl, $6 = — , H , (A.2) 

and obey the hrst law of black hole thermodynamics 

dM = TbdSb + %dQ + VbdP. (A.3) 

For sufficiently small pressure 0 < P < Pc = l/967rQ^, the Gibbs free energy 

G = M- TbSb. (A.4) 
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Gi 


T2 


Ti 


Figure 11. Behavior of Gi. The behavior of the inner horizon Gibbs free energy Gi is 
displayed for P = 0.003 and Q = 1. The two black curves correspond to thermodynamically 
stable inner horizon configurations whereas the dashed red curve in between Ti and T 2 cor¬ 
responds to the thermodynamically unstable configuration of the inner horizon for which the 
outer horizon maintains temperature Tq and undergoes a first order phase transition. 


possesses a swallow tail, which represents the gravitational analogue of a liguid/gas 
phase diagram, see Fig. 10. Namely, there is a hrst order small black hole/large black 
hole phase transition. During this transition, the radius of the black hole suffers 
from a jump while the temperature Tf, and pressure P remain constant. 

However, the black hole spacetime also possesses an inner black hole horizon, lo¬ 
cated at Tj, 0 < Tj < Th, f{ri) = 0. What happens to this horizon during the above 
phase transition? It is straightforward to show that the inner horizon can be assigned 
the following thermodynamic quantities: 


Ti = 


Airvi 


1 -b SirPrf — 




Si = vrr; 


= 


Q 

Ti 


Vi = 


4 3 

-7rr, , 

3 * 


(A.5) 


that formally satisfy the following hrst law: 


dM = -TidSi + ^idQ + VidP , 


(A.6) 


Identifying the enthalpy Hi = —M we can write this as dHi = T^dSi — VidP — ^idQ 
(note the inevitable minus sign in the work terms) and hence identify the corresponding 
Gibbs free energy as 

G, = -M- T,S, . (A.7) 

The behavior of Gi is displayed in Fig. 11. We observe a monotonic behavior seemingly 
without any phase transitions. This is misleading however. It is easy to show that 
while the black hole undergoes a hrst order phase transition and the horizon radius 
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Figure 12. Thermodynamic description of both horizons. The Gibbs free energy 
G is displayed for both positive and negative temperatures to capture the thermodynamic 
behavior of the inner and outer horizons simultaneously. We have set Q = 1 and P = 0.001; 
for such a choice T 2 becomes too large to be displayed in this hgure. 

jumps, the inner horizon jumps as well. Consequently, the inner horizon temperature 
Tj (which is a monotonic function of rj) also jumps from Ti to T 2 , unlike the black hole 
horizon temperature T;, = Tq which remains constant. This means that in fact only a 
part of the displayed Gibbs free energy Gi corresponds to a stable equilibrium state 
and the region in between Ti and T 2 is thermodynamically unstable. In other words, 
the black hole hrst order phase transition leaves an imprint on the inner horizon: the 
inner horizon suffers from both a jump in radius and a jump in temperature. 

Let us hnally remark that we can write 

Gi{ri,P,Q) =ri,P,Q), Ti{ri, P,Q) =-Tb{n-i-Vi, P,Q). (A.8) 

For this reason we may formally obtain the thermodynamic information about both 
horizons by studying G, (A.4), extended to all admissible non-negative radii and hence 
formally to negative temperatures. Negative temperature corresponds to the positive 
inner horizon temperature, while in such a region the thermodynamic equilibrium state 
maximizes (rather than minimizes) G. This is displayed in Fig. 12. 

B Rotating de Sitter black holes and their thermodynamics 

In this appendix we review d-dimensional Kerr-de Sitter black holes and their thermody¬ 
namics [20, 43-45]. The metric solves the vacuum Einstein equations with cosmological 
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constant A, (1.4), and reads 


9m / , 
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. + af r" , ^ 

1=1 ‘ k=l 


% 

P 


(B,2) 


Here TV = [(d — l)/2], where [A] means the integer part of A, and we have dehned e to 
be 1 for even d and 0 for odd d. The coordinates /Xj are not independent, but obey the 
constraint 


N 


+ = 1 - 


(B.3) 


k=l 


We have up to three horizons (two black hole horizons and a cosmological dS horizon). 
These occur at positive radii determined by /(r) = 0. 

B.l Even dimensions 

In even dimensions {d = 2N + 2), we hnd the following thermodynamic quantities 
[20, 45]. The ‘mass’ and angular momenta read 


M = 


mud-2 


E^. J- 


maiUd-2 


where Ud -2 stands for the volume of the unit-radius (d — 2)-sphere, 


(B.4) 
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The cosmological horizon is assigned the following entropy, temperature, angular ve¬ 
locities, and thermodynamic volume: 
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Similarly, for the black hole horizon and the inner horizon we have the following quan¬ 
tities: 
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B.2 Odd dimensions 


In odd dimensions (d = 2N -|- 1) the mass, temperature, and entropies get modihed as 
follows: 
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All other quantities (fl’s, J’s, and Id’s) remain as in the even-dimensional case. It is 
easy to verify that in both even and odd dimensions these quantities satisfy the hrst 
laws (1.1)-(1.3) as well as the Smarr formulae (1.7)-(1.9). 













































The asymptotically AdS black holes are formally obtained by setting I ^ il. The 
corresponding thermodynamics in extended phase space has been studied e.g. in [8-11]. 
In particular it has been shown that for d < 5 one only observes the Van der Waals- 
like phase transition, similar to the spherically symmetric charged AdS case. More 
interesting is the behavior of black holes in d > 6 dimensions. In particular, in d = 6 
dimensions one can observe the reentrant phase transitions [8], triple points [9], or Van 
der Waals behavior, dependent on the ratio of the two angular momenta q = Ji/J 2 - 
Interestingly, as we have seen in the main text, the reentrant phase transition for d = 6 
dimensional black holes can also occur for the dS case. 
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